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The experimental data of the COSY-GEM collaboration for the recoil- 
free transfer reaction p(^^Al,'^He)7r~p'X, leading to the formation of bound 
state of eta (ry) meson in ^^Mg nucleus, is reanalyzed in this paper. In 
particular, predicted values of binding energy and half- width of the 77-mesic 
nucleus ^^Mg,,, given by different theoretical approaches, are compared 
with the ones obtained from the experimental missing mass spectrum. It 
is found that the spectrum can be explained reasonably well if interference 
effect of another process, where 77 is not bound in ^^Mg but is scattered 
by the nucleus and emerge as a pion, is taken into account. The data also 
indicate that the interaction between N*(1535) and a nucleus is attractive 
in nature. 

PACS Numbers: 24.10.Eq, 24.10.-I, 21.10.Dr 

1. Introduction 

The existence of eta-mesic nucleus, a bound state of r] meson in a nuclear 
orbital, was first predicted by us in 1986 [1]. The formation of such a nucleus 
is a consequence of the attractive //-nucleon (r/N) interaction in the threshold 
region [2]. This attractive nature of the rjH interaction is due to the fact 
that the N*(1535) resonance is situated not too far above the threshold of 
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the r/N channel (~ 1488 MeV). Our prediction has been corroborated by 
other researchers and can be found in the vast amount of theoretical work 
published during the past 20 years [3, 4J. 

Experimental efforts to detect r/-mesic nucleus, on the other hand, did 
not yield positive results [21 E]- It was only recently the COSY-GEM col- 
laboration [7] confirmed unambiguously the detection of mesic-nucleus in 
the recoil-free transfer reaction p(^^Al,^He)7r~yX. The experimental kine- 
matics was chosen such that in the first stage of the reaction the rj was 
produced nearly at rest, thus favoring its capture by the residual nucleus 
^^Mg to form the ry-mesic nucleus that we denote as ^^Mg^: 

p + 27A1 -^{rj + 25Mg) + 3He ^ ^s^g^ ^ 3jjg_ 

Because of energy conservation, the bound r] cannot reappear as an ob- 
servable in the decay products. Rather it interacts with a target nucleon 
resulting in the emission of a nearly "back-to-back" 7r~p pair in the labora- 
tory, i.e., 

25Mg^ (r? + n) + 24Mg (vr" + p) + X. 

We may call the above multi-step reaction Process M as it proceeds through 
the doorway state containing mesic nucleus. By fitting the missing-mass 
spectrum with + where and fg are, respectively, the back- 

ground amplitude and a Gaussian function, the COSY-GEM collabora- 
tion has determined the values of the binding energy and FWHM of the 
mesic nucleus ^^Mg^ to be, respectively, e = (—13.13 ± 1.64) MeV and 
r = (10.22 ± 2.98) MeV (or r/2 ~ 5.1 ± 1.5 MeV. We refer the readers to 
ref. [7] for the experimental details and the accompanying analysis of the 
data. 

The experimental verification of the existence of mesic-nucleus by COSY- 
GEM collaboration has opened up new avenues in the study of exotic nuclei. 
In this work, we reanalyze the data and investigate how theory measures up 
to the experimental results. In particular, we want to address the following 
questions: 

1. Are the published values of s-wave ??N scattering length {urjN) consis- 
tent with the data? 

2. Can the data be described within the context of the models and for- 
malism used by us [8| l9], which predicted the existence of r/-mesic 
nucleus? 

3. Is it possible that under the conditions the experiment was done, r/ 
got captured in the excited state of ^^Mg? 
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4. Is there any other process that could have contributed to the observed 
experimental spectrum? 

5. Can the observed spectrum provide information about the nature of 
the N*(1535)-nucleus interaction? 

We believe that answer to the above questions will enhance our understand- 
ing of ?7-nucleus reaction dynamics and properties of ry-mesic nucleus. 



2. Theoretical analysis 

To obtain the binding energy and half-width of r/-mesic nucleus, we solve 
the relativistic three-dimensional (covariant) integral equation 

— V'(k') + j dk < k' I y I k > V(k) = Ei){\s!) , (1) 

where k and k' are the momenta of rj in the jy-nucleus cm. frame, and /x 
is the reduced mass of the ?7-nucleus system. The eigenvalue E is complex 
and is given by ^ = e — iT/2 with e < and F > 0. 



2.1. Scattering-length approach 

In the scattering-length approach, the first-order low-energy ry-nucleus 
optical potential is 

where /(k — k) is the nuclear form factor, iW^ is the mass of r] meson, 
and Mjv is the nucleon mass. The only input in calculating e and T /2 is, 
therefore, the scattering length (i^^n- This approach can be also termed as 
"on-shell" calculation because the scattering length is defined at jyN on-shell 
threshold energy -|- Mjv. 

The value of a^jv that reproduces the e and r/2 determined by the 
COSY-GEM analysis equals to (0.292+0.077i) fm. We have noted, however, 
that all published theoretical models give a^jv having an imaginary part 
that is at least a factor of 2.4 greater than 0.077 fm. The existence of 
this significant difference is theoretical-model independent and needs to be 
understood. In the next subsection we analyze the data more thoroughly by 
using a microscopic optical potential which takes into account the energy 
dependence of 77N interaction below the 77N threshold. 
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2.2. Microscopic (off-shell) optical potential approach 

The momentum-space matrix elements of the microscopic r/-nucleus op- 
tical potential is [H [8] 

< k'|F|k >=< K'|t^7v(VF)|K > . (3) 

The variables k, and k' are the initial and final ryN relative momenta. The 
total energy of the rjH system in its center-of-mass frame is 

W = M^ + Mn+<Bn>, (4) 

where < > (< 0) is the average binding energy of the nucleon. We have 
calculated the matrix elements of the scattering operator tr,NiW) by using 
the coupled-channel isobar (CCI) model of Bhalerao and Liu [2], which has 
its parameters determined from fitting the ttN Sh phase shifts. The same 
CCI model was used in all of our earlier calculations that predicted the 
existence of r/-mesic nucleus. More importantly, the model provides us with 
the detailed energy dependence of the interaction in the energy region where 
r] can be bound in a nucleus. It contains strong-interaction form factors and 
satisfies off-shell unitarity. 

Because for bound-state problems the contributions of the p- and d-wave 
r/N interactions are negligibly small, only the s-wave r/N interaction needs 
to be taken into account. The t-matrix of the CCI model, therefore, has the 
form 

< K'\tr^N{W)\K >= Kv{K',A)AiW)v{K,A), (5) 

where K is a. kinematic factor, v is the s-wave off-shell form factor, and A 
is the range parameter. The energy-dependent amplitude is 

where g is the r/NN* coupling constant. The expression for D{W) is 

D{W) = - {M° + VN*iW) + ^f''^^{W)+iIm[^''^'iW)]}, (7) 

where is the bare mass of the resonance, and the N*(1535)-nucleus 
interaction is defined as 

Vn'{W) = r{W)+ne[^''^%W)]. (8) 

In the above expressions, ^•^''^^(VF) is the complex self-energy of N* arising 
from its decay to the ?7N, vrN, and vrvrN channels. The self-energy arising 
from absorption or annihilation of the rj and vr resulting from N* decaying 
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into r?N, vrN, and vrvrN is denoted by Tf'^''{W). The r{W) in eq.® is a real 
quantity which accounts for ah other subthreshold N*-nucleus interactions. 
Although has been evaluated within the framework of local 

density approximation [TO], theoretical models for r{W) as well as for cal- 
culating 7^e[S'^^''(^^)] that can be systematically checked against data are 
still unavailable. We, therefore, treat Vn*{W) as a parameter. 

We emphasize that unitarity of an optical potential requires that the 
imaginary parts of Tj°'^'^{W) and T,-!'^'^'^ (W) should have the same sign. Con- 
sequently, I Im[D{W)] \>\ Xm\Ti^'^'^^{W)\ \. In our calculations, careful 
attention is paid to this requirement. 

Upon applying the CCI model to ^^Mg with Vn*{W)=^ and < Bjq >= 
—30 MeV, we obtained binding energy e = —6.5 MeV and half-width r/2 = 
7.1 MeV. (We point out that < Bn >= -30 MeV is based on the findings 
obtained from extensive pion- nucleus studies in the literature.) As we can 
see, although the calculated half-width is comparable with the data, the 
magnitude of the calculated binding energy is about 7 MeV too small. 

To reproduce the COSY-GEM values of e = -13 MeV and r/2 = 
5.0 MeV requires that the amplitude A = -(0.0521 + 0.0099i) fm^. A 
knowledge of the value of A allows us to determine W ^ and hence < Bjsf >, 
by means of two equations for Im[D{W)] obtained from eqs.® and dTj). 
They are 

Mowl = (|,)x™(^). P) 

and 

Im[D{W)] = -{Im[i:f'^%W)] + Im[i:'''" (W)]} . (10) 

With the fitted value of A and the calculated imaginary part of the self 
energies, we find that the values Im[D{W)] given by the last two equa- 
tions are same if = 1125 MeV. This implies < -Btv >= —360 MeV, which 
is clearly unrealistic. We conclude, therefore, that a simple application of 
optical potential does not yield satisfactory result, and a reexamination of 
the whole reaction is in order. 

2.3. Excited state of^^Mg 

In view of the above findings, it is reasonable to investigate whether 
the seemingly large binding energy of ^^Mg,, could be due to the binding 
of T] on excited states of ^^Mg. To answer this question, we recall that 
the experimental spectrum as a function of the binding energy E oi r] is 
deduced from measurements of the missing mass AM. In terms of the 
mass M25 and excitation energy Ex (> 0) of ^^Mg, energy conservation and 
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recoil-free kinematics give 

AM = M25 + Mr, + E^ + E. (11) 

With known values of the masses, we find that the experimental centroid of 
the spectrum at AM = 23.803 GeV [7] leads to 

E = -{13 + E^) (MeV). (12) 

From this equation, we see that the centroid of the observed spectrum at 
— 13 MeV implies = 0. We, therefore, conclude that the r/ was captured 
in the ground-state of ^^Mg and that the noted difference between theory 
and experimental data cannot be due to r] being bound on excited states. 



3. Two doorway processes 

The inability of the theoretical models to explain the data within the 
framework of doorway process M alone leads us to investigate whether there 
is another process contributing to the observed spectrum but has been over- 
looked in the above discussions. Indeed, the r] produced in the intermediate 
state can scatter from the residual nucleus and emerge as a pion, without go- 
ing through formation of an r/-mesic nucleus. This latter multi-step process 
is 

p + 2'^Al ^ (r? + ^^Mg) + ^He 



followed by 



■q + ^^Mg {r] + n) + ^^Mg (vr" + p) + X. 



We call this Process S (S for scattering). It should be emphasized that 
because the Processes M and S lead to the same final experimental state, 
the effect of interference between them must be taken into account. We, 
therefore, fit the experimental spectrum with the function a\fs + /mP, 
where fs and fu are the amplitudes for the Processes S and M, respectively, 
and a adjusts the overall magnitude of the spectrum. The amplitudes are 
parameterized as 

The parameters A and 9 represent the relative strength and phase between 
the two processes at £" = e. The energy dependence of fs is neglected in 
the present phenomenological approach because it is much smoother than 
/a/. The half-width controls the sharpness of the peak's structure of the 
calculated spectrum. With the aid of eq. (|13|) . we carried out fits to the 
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Fig. 1. Spectra obtained with e - ir/2 fixed at : (a) -(6.5 + 7.1i) MeV (dotted 
curve); (b) -(8.0 + 9. 6i) MeV (dashed curve); (c) -(10.0 + 14.0i) MeV (sohd curve). 
The data are from ref. [7j. 



Table 1. Fitted values of the parameters. 



Fit 


e 


r/2 


a 


A 


e 


w 


<Bn> 


Vn* 


(a) 


-6.5 


7.1 


1.83 


1.75 


0.62 


1458 


-30 





(b) 


-8.0 


9.6 


1.91 


1.14 


0.79 


1458 


-30 


-15 


(c) 


-10.0 


14.1 


2.54 


1.17 


1.06 


1458 


-30 


-42 



(e, r/2, W, < Bn >, and Vn* are in MeV, a is in counts/4-MeV, A is 
dimensionless, and 9 is in radians.) 



data by varying a, A, and 6 while keeping fixed the values of e and r/2, 
which were calculated using the microscopic optical potential. The fits for 
three different sets of parameters are shown in figure 1. The values of the 
parameters are listed in tabled) 

A few remarks about the fits are in order. First, fit (a) is obtained 
without nuclear medium effect (i.e., Vn* = S'^''*(VF) = 0, as in our previous 
predictive calculations), but with the interference effect between fs and 
/m included. As can be seen, interference shifted the observed peak from 
—6.5 MeV to —9.0 MeV. Fits (b) and (c), on the other hand, are obtained 
with Vj\f* < 0. One can see that the final positions of the peak of q|/s + /a./P 
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are approximately at —11 MeV for fit (b) and —13 MeV for fit (c). In 
general, the interference effect causes the final peak positions, as indicated 
by the respective curves in figure 1, to move closer to the experimental 
centroid. In fact, in the case of fit (c) it is equal to the experimental centroid. 
The shift in the position of the peaks clearly indicates the importance of 
the effects of interference. The negative value of V/v* for fits (b) and (c) 
indicates that the data favor an attractive N*(1535)-nucleus interaction. 

One may ask why we did not fit the data by using larger |e| and smaller 
r/2 or stronger Viv*? Without going into details, we just mention that these 
quantities cannot be treated as completely free parameters. As pointed 
out in section 12.21 their relative strengths are constrained by the unitarity 
requirement of the optical-potential: | Im[D{W)] \> \ Zm[S'^'"'^^(VF)] |. 



4. Conclusion 

Our analysis shows that two reaction processes are contributing to the 
observed spectrum of bound t] in ^^Mg. The quantum interference between 
the two processes results in a binding of r] that is weaker than suggested 
by the centroid of the observed spectrum. The present analysis gives the 
binding energy between —8 and —10 MeV and the half- width between 10 
and 14 MeV if ^^Mg is in its ground state. 

The present analysis also indicates that the real part of interaction be- 
tween N*(1535) and a medium-mass nucleus is attractive at energies below 
the r/N threshold. This latter information should be of value to nuclear 
physics studies involving the baryon resonance N*(1535). 

Microscopic calculation of fs and Jm is in progress, where we calculate 
both amplitudes by using the same r/- nucleus optical potential. Accordingly, 
there is no need of the phenomenological parameters A and 9. Preliminary 
results reaffirm the basic findings reported here, namely, (i) the existence 
of two reaction processes, (ii) the interference effect between them causes 
the peak oi \ fs + fu P to appear at a binding energy much stronger than 
the actual binding energy of ^^Mg^, and (iii) the attractive nature of the 
N*-nucleus interaction. Details will be reported in a future publication. 

One of us (Q.H.) would like to thank Dr. Pawel Moskal for the hospi- 
tality extended while he was in Krakow, Poland to attend the International 
Symposium on Mesic Nuclei. 
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